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Abstract 

A model for highly non-linear transport in heterogeneous media consisting of anisotropic particles 
with a preferred orientation is proposed and applied to the current transport in weakly textured 
magnesium diboride, MgB2- It essentially explains why, unlike in conventional superconductors, 
a significant macroscopic anisotropy of the critical currents can be induced by the preparation of 
MgE>2 tapes. The field and angular dependence of the critical current is calculated for various 
degrees of texture and compared to experimental data. 



PACS numbers: 74.25. Sv,74.81.Bd,74.70.Ad,81.40.Ef 



I. INTRODUCTION 



The critical current in polycrystalline MgB2 is a particular example of nonlinear transport 
in heterogeneous media 1 . The intrinsic upper critical field anisotropy results in different 
critical current densities in differently oriented grains, when a magnetic field is applied 2 . 
The problem was solved by integration of the (local) current densities over the percolation 
cross section, for which the distribution function of the local critical currents has to be 
known. Although the local currents can be described by standard pinning models 3 -, the 
resulting macroscopic currents are significantly altered. The field dependence of the critical 
current increases and the current becomes zero well below the upper critical field 2 . 

In the present letter, a heterogeneous and macroscopically anisotropic system is con- 
sidered, i.e. the distribution function of the relevant transport property (e.g., the critical 
current density J c ) depends on its orientation. This is realized in weakly (or partially) tex- 
tured MgB2 and explains the Jc-anisotropy^&^&^ii^, which can exceed the intrinsic 
anisotropy of MgB 2 (5-6) 13 by orders of magnitude. The J c -anisotropy can be induced by 
some preparation techniques of MgB 2 tapes, which was impossible in the conventional super- 
conductors NbTi and NbaSn since they are intrinsically isotropic. On the other hand, high 
temperature superconductors must be highly textured for enabling large transport currents. 
Thus, MgB 2 is the first superconductor needing a model for partial texture as description. 
The distribution function of J c in MgB 2 will be derived in the following, but it is straightfor- 
ward to apply the model to any other highly nonlinear transport phenomena in anisotropic 
heterogeneous media. 

II. EXPERIMENTAL 

Powders of magnesium (Goodfellow, 99.8%) and amorphous boron (Starck, 95-97%) 
were milled together for 50 hours at 250 rpm in a planetary ball mill (Retsch PM 4000) 
under protecting argon atmosphere. The powder was filled into iron tubes and deformed by 
swaging and rolling into tapes and heat treated. One tape was doped with nanostructured 
carbon (99.999%), which had been added to the precursor powders before milling (nominal 
MgBi.g7C0.13). 

The critical current was measured by direct transport measurements in liquid helium. 
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FIG. 1: Spherical coordinate system. The magnetic field is assumed to be parallel to the z-axis. 9 
and ip denote the polar and azimuthal angle, respectively. The preferred orientation (texture) of 
the grains lies in their c-axis parallel to t. 

The texture was characterized by measuring rocking curves of the MgE>2(002) reflection using 
high energy synchrotron X-ray (E = 77.5 keV) diffraction. Details of all measurements and 
the sample preparation are given elsewhere^, since we will focus on the proposed model, 
while the experimental data are used for comparison. 

III. MATHEMATICAL MODEL 

The macroscopic critical current density in anisotropic polycrystalline superconductors 
(or any other transport property in heterogeneous media, which hardly depends on the 
driving force) can be calculated from 2 



The percolation threshold p c is defined as the minimum fraction of superconducting grains 
for a continuous current path. p(J) is the fraction of grains with a (local) J c exceeding J 
and J max is given by the condition p( J max ) = p c . The field dependence of the (local) critical 
currents is assumed to follow the well established model for grain boundary pinning 3 




dJ. 



(1) 




(2) 
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for B < B C 2 and J ° = otherwise. For the sake of simplicity, any difference between 
the upper critical field, B c2 , and the irreversibility field, B irr , is neglected. The angular 
dependence of the critical current is modeled by the scaling approach proposed by Blatter 
et al.— , i.e. B is multiplied by J cos 2 9 + 7 -2 sin 2 6, where 7 denotes the anisotropy of the 
upper critical field B^/B^ and 9 the angle between the applied magnetic field and the 
crystallographic c-axis; thus, B c2 in Eq. (j5J) corresponds to B° 2 , the upper critical field for 
B\\c. 

The angular distribution of the grain orientation, p a (9), is needed for calculating p(J). 
Each grain is represented by a point on the surface of the unit sphere (r in Fig. [TJ, which 
corresponds to the orientation of its c-axis. The density of these points, / a (0, <p) is per 
definition constant for a random grain orientation. The number of grains with a polar angle 
equal to or above 9 is obtained from integration of f a (9, p) over the corresponding surface 
area 

Pa(0) = I" r W, V )dy sin0W, (3) 
Je Jo 

which reduces to 27r i / a cos# for constant / a . The sine of 9' is the radius of the circle, along 
which the integration over p> is carried out (cf. Fig. [1]). Due to the symmetry of the problem, 
the integration can be restricted to the upper half of the unit sphere. The relative fraction 
of such oriented grains is obtained by setting / a = 1/27T (the surface of the unit sphere is 
4n). This leads to p a (9) = cos9 and to the correct result p a (0) = 1 for a random grain 
orientation. 

Partial texture can be introduced by a varying density f a (9, <p). For uniaxial texture, the 
density becomes independent of the azimuthal angle (ip), if the preferred orientation of the 
grains' c-axis is assumed to be parallel to the z-axis. A natural choice for the variation along 
the polar angle is a Gaussian (or normal) distribution 

q2 

f a {9,p)ccexp{-^). (4) 

with standard deviation ct t . The proper normalization is obtained as above by integration 
of fa,{9,ip) over half the unit sphere and p & {9) can then be calculated from Eq. The 
chosen density function was confirmed by synchrotron x-ray diffraction with a t = 29.3°—, 
but it is straightforward to use any other distribution. 

If the direction of the applied magnetic field is chosen as the z-axis, the preferred ori- 
entation of the grains' c-axis is in general not parallel to this direction. The situation is 
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illustrated in Fig. [TJ The preferred orientation is defined by t with spherical coordinates 
(6> t , <p = 0). The angle, a, between t and an arbitrary direction r is given by 

cos a = abs(sin 9 cos tp sin 9 t + cos 6 cos 9 t ) . (5) 

The density function / a (#, <p) is obtained from Eq.(Tj| by replacing 9 by a. It depends on 
the polar and azimuthal angle and the two dimensional integration ([3]) has to be calculated 
numerically 

At fixed magnetic field, p(J) can be calculated from p&(9), since J,l ocal depends monoton- 
ically on 9. Finally, J™ acro is obtained by numerical evaluation of integral ([1]). 

IV. RESULTS AND DISCUSSION 

The relative abundance of grains, whose c-axis is tilted by an angle 9 from the direction of 
the texture t, is obtained by integration of / a (#, V 9 ) over V 9 ; which is the negative derivative of 
Pa,(9). It is plotted in Fig. [2] for various angles a t , which characterize the degree of texture. 
With increasing texture, the maximum shifts from n/2 for randomly oriented grains to 
smaller angles. (In the limit a t —>■ (perfect texture), the distribution function becomes a 
5-function.) For weak texture, the majority of grains is not oriented with the c-axis close to 
the direction of the uniaxial texture, but the fraction of grains with perpendicular orientation 
is reduced and becomes negligible for a t ~ 25°. 

The influence of partial texture on the field dependence of the critical currents is shown 
in Fig. [3j The two panels refer to different anisotropies; 5 is typical for clean MgB2, 3 for 
moderately dirty MgB 2 .— The magnetic field is normalized by B°%. B c c2 is I/7 (0.2 and 0.3). 
The percolation threshold was fixed to 0.25, which is an average value in polycrystalline 
MgBz^^iSiil. The central bold line was calculated for the untextured material, the graphs 
below and above this line refer to the parallel (9 t = 0°) and perpendicular (9 t = 90°) 
field orientation. With decreasing a t (increasing texture) the critical current becomes more 
anisotropic and the data converge much faster to the limiting case of perfect texture (thin 
solid lines) for the parallel orientation, where the difference between perfect and partial 
texture becomes small at a t ~ 20°. The J c -anisotropy decreases with decreasing 7 but the 
qualitative behavior is the same in both panels. 

The angular dependence of the critical currents is plotted in Fig. HI The assumed texture 
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FIG. 2: (Color online) Grain orientation with respect to the direction of partial texture. The 
maximum shifts continuously from the perpendicular orientation (90°) for random grain orientation 
to zero in the limit of perfect texture (at — * 0). 

is stronger in the lower panel, the anisotropy (7 = 5) and p c (0.25) were kept constant. At 
low magnetic fields, the J c -anisotropy is small for the more weakly textured material (1.2 at 
B/B C 2 = 0.01), but increases with texture (to yfy for B — > for perfect texture, line graph 
in Fig. H]). The J c -anisotropy monotonically increases with magnetic field until it diverges 
for any a t < 00, when the currents reduce to zero in the parallel orientation (6 t = 0°). This 
was observed experimentally 5 - 1 ^. (Note that t is perpendicular to the surface of MgB 2 tapes, 
hence the parallel orientation (9 t = 0°) corresponds to the field perpendicular to the tape.) 
At fixed 7, pc, and a t , the J c - anisotropy is only a function of the reduced magnetic field 
b = B/B c2 , which explains its temperature dependence at fixed magnetic field: b increases 
with temperature, since B c2 decreases. 

Experimental data (symbols) for the critical current are compared to the prediction of the 
present model (line graphs) in Figs. [5] and EJ The agreement is excellent, given the numerous 
assumptions/simplifications (pure grain boundary pinning^, anisotropic scaling 15 , Gaussian 
distribution etc.). The correct choice of the (unknown) parameters is difficult, since a 5 
parameter fit (B c2 , 7, a t , p c , and a constant prefactor, J c0 ) does not lead to significant 
results. Therefore, p c was fixed to 0.25 and a t to 29.3°—. The remaining parameters were 
adapted to fit the experimental J C {B) data for both field orientations simultaneously. B c2 = 
23.2 (20.5) T, 7 = 4.05(1.8) and J c0 = 2840 (860) A were chosen for the undoped (doped) 
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FIG. 3: (Colour online) Influence of partial texture on the critical currents. The given angles refer 
to at- The central bold lines represent untextured MgB2- The dotted lines below and above the 
bold lines correspond to the field orientation parallel and perpendicular to the direction of texture. 
The limiting curves represent perfect texture. 

tape. The parameters of the doped tape (T c = 29.5 K) agree fairly well with results on a 
carbon doped (9.5%) single crystal with a similar T c = 30.6 K 18 , where 7 slightly below 2 
and B C 2(0) = 25 T (suggesting a value of around 21 T at 4.2 K) were reported. The higher 
anisotropy of the undoped tape is expected from its higher T c (34.4 K) and also B c2 tends 
to decrease for transition temperatures of below ~ 32 K 13 . The larger I cQ of the undoped 
tape confirms that carbon doping decreases pinning^ 1 ^ and the connectivity. 

A discrepancy between experimental and calculated data was found at small currents 
(< 1 A, 6 t = 0°, not shown in Fig. [5]), which might be a consequence of the finite voltage 
during transport measurements (1/iVcm -1 ), which allows small currents passing normal 
conducting grains. Such normal currents are excluded in the present model and might add 
a significant contribution at small critical currents (by connecting otherwise disconnected 
superconducting clusters), but they can certainly be neglected at high critical currents. On 
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FIG. 4: (Color online) Angular dependence of the critical current. 7 = 5 and p c = 0.25 are 
assumed. The texture is stronger in the lower panel. The graphs without symbols refer to perfect 
texture (anisotropic scaling approach 15 ) at B/B C 2 = 0.01. The critical current anisotropy increases 
with field. 
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FIG. 5: (Color online) The model (line graphs) successfully describes the measured critical currents 
(symbols) for both main field orientations. 
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FIG. 6: (Color online) Calculated (lines) and experimental (symbols) angular dependence of the 
critical currents at 7.5 T in liquid helium. 

the other hand, the discrepancy could also result from inhomogeneities, texture gradients^, 
or a different distribution function of the grain orientation. Note that the Gaussian distri- 
bution function was also used in Ref. (] and an even higher degree of texture (a t ~ 20° and 
~ 25° according to our definition) was found near the sheath of cold rolled ex-situ tapes, 
decreasing towards the center. 

The angular dependence of J c was calculated with the same parameters, which describe 
the data very well. However, without a knowledge of the anisotropy the texture (a t ) cannot 
be derived due to the mutual influence of a t and 7. 

V. CONCLUSIONS 

The proposed model successfully describes the field and angular dependence of the crit- 
ical currents in weakly textured MgB 2 tapes and explains the large difference between the 
intrinsic and extrinsic (J c ) anisotropy. 
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